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ABSTRACT 

We present an estimate of the mass of the supermassive black hole (SMBH) in the nearby type-1 
Seyfert galaxy NGC 1097 using Atacamma Large Millimeter/Submillimeter Array (ALMA) observa¬ 
tions of dense gas kinematics. Dense molecular gas dynamics are traced with HCN(J = 1 — 0) and 
HCO + ( J = 1 — 0) emission lines. Assuming a host galaxy inclination of 46°, we derive a SMBH mass, 
Mbh = 1.40t°;£ x 1O 8 M 0 , and an I-band mass to light ratio to be 5.141 q’ 04 , using HCN(J = 1-0). 
The estimated parameters are consistent between the two emission lines. The measured SMBH mass is 
in good agreement with the SMBH mass and bulge velocity dispersion relationship. Our result show¬ 
cases ALMA’s potential for deriving accurate SMBH masses, especially for nearby late-type galaxies. 
Larger samples and accurate SMBH masses will further elucidate the relationship between the black 
hole (BH) and host galaxy properties and constrain the coevolutionary growth of galaxies and BHs. 
Subject headings: galaxies: supermassive black holes - galaxies: individual (NGC 1097) - black hole 
physics 


1. INTRODUCTION 

Recent observations suggest that supermassive black 
holes (SMBH) resid e in the centers of most mas¬ 
sive galaxies (e.g., iKormendv fe Hoi I2013L and refer¬ 
ences therein). In the nearby universe a variety of 
host galaxy properties are known to be correlated 
with the central SMBH mass. For instance there 
is a tight correlation between the SMBH mass and 
the bulge luminosity (Mrh — Lhnige relation , e.g., 
Kormendv fe R ichstone 1995]; IMcConnell fe Mai 120131 : 
Kormendv fe Ho 201 3]), the bulge ma s s (Mbh — Mhnip-e 
relation, e.g.. iM agorrian et al.lfl9 98: Marconi & Hunt 
120031: iBeifiori et al. 120121: IKormendv fe Hoi l2013h , and 
the central velocity dispersion (Mrh — cr re l ation, 
e.g.. Ihbr rarese & Merritt 120001: Giilteki n et al.l 120091 : 
IMcConnell fe Mai 120131: IKormendv fe Hoi[20131 ). These 
empirical correlations suggest that black holes (BHs) 
may play a key role in the growth and evolution of galax¬ 
ies. 

Recent studies, however, reveal that the correlation 
between the SMBH mass and bulge/host galaxy prop¬ 
erties ar e more complex than o riginally thought. For 
instance IMcConnell fe Mai (J2013D showed that the best- 
fit Mrh — cr relation differs between early- and late-type 
galaxies. The difference results in the SMBH mass to 
be two times larger at a given velocity dispersion for the 
early-type Mbh ~ cr relation than for the late-type sys¬ 
tems. In contrast, there does not seem to be a correlation 
betwee n galaxies hosting psued o-bulges and their SMBH 
masses (jKormendv fe Holl2013[ ) . Much of the uncertainty 
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comes from the scarce number of measurements of the 
SMBH masses especially in late-type galaxies, which is 
especially difficult with current methods for measuring 
the SMBH mass. 

The most reliable way to estimate SMBH mass is to 
use dynamics/kinematics of gas or stars near the SMBH 
in a galaxy. So far, SMBH masses have been mea¬ 
sured in the following wa ys: using proper motion of stars 
around the SMBH (e.g., iGenzel et al.l 1200(1 iChez et all 
mm using proper motion and dynamics of megamasers 
(e.g., iMivoshi et al.l 119951: iKuo et all l2011f ). stellar dy¬ 
namics (u sually only for elliptical or lenticular galax¬ 
ies) ( e.g.. [Dr essier & Richston^ 119881 : IMcConnell et al.l 
20111 lOnken et al.ll20Tl. and ionized gas dynam ics (e.g., 
Macchetto et al.l 119971: I de Francesco et al.ll2008l ). 

The method using the proper motion of stars around 
the central SMBH is the simplest of the four, but spa¬ 
tially resolving stars around a SMBH, besides the one 
in our Galaxy, is not currently possible. The SMBH 
mass measurements in galaxies using megamasers, which 
are rare and difficult to find, requires observations with 
very high angular resolu tion (~ 0.3 milliarcsecond, for 
example, in Mivos hTet al.l Q.995), accomplished with 
Very Long Baseline Interferometer. Stellar dynamics 
have been used to measure the largest number of SMBH 
masses at this point. An orbit superposition model 
(Schwarzschild 1979j) is fit to spectroscopic observations 
and the method is primarily restricted to elliptical and 
lenticular systems. The errors in this method can be 
amplified when a mass profile is asymmetric. Ionized 
gas dynamics can be observed in larger samples (at least 
compared to the use of stellar or maser dynamics), by 
taking spectra of a galaxy using multiple slits or Inte¬ 
gral Field Units (IFUs). The weak point of this method 
is that ionized gas is not necessarily settled into a pure 
rotating disk because it is more easily affected by non¬ 
circular motion from turbulence, shocks, radiation pres¬ 
sure, outflows, etc. Inferring the velocity field from slit 
observations is also problematic because such measure¬ 
ments may not always show non-circular motions which 
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could be affecting the gas dynamics. Note that recent In¬ 
tegral Field Unit (IFU) observations can avoid this prob¬ 
lem and are improving to be very useful to obtain velocity 
fields. 

Deriving SMBH mass from mo lecula r gas dynamics 
was not accomplished until iDavis et~ah (1201 3) because 
millimeter-wavelength interferometers did not have suf¬ 
ficient angular resolution and sensitivity to map out the 
precise kinematics of molecular gas around a black hole. 
The SMBH mass in the central region of NGC 4526 
was measured using the observed CO emi ssio n line from 
its circumnuclear, molecular-gas dis k (IDavis et all 2013b 
This method was also proposed bv lOnishi et all ( 2013 ). 
This technique is similar to the one used for ionized 
gas dynamics except that a molecular gas velocity field, 
which is primarily observed at millimeter/submillimeter 
wavelengths with interferometers, can provide very high 
spectral resolution with full two-dimensional kinematics 
compared to optical or near-IR observations. Another 
major strength for molecular gas observations is that 
they can be utilized across the Hubble sequence from 
early- to late-type galaxies, as long as there is molecu¬ 
lar gas in the galaxy centers. Thus molecular gas ob¬ 
servations with new facilities like Atacamma Large Mil¬ 
limeter/submillimeter Array (ALMA) offer a new and 
promising avenue for increasing SMBH measurements^ A 
figure of merit for this method is summarized by iDavisI 

dmTl . 

In this pam per, we use the method initially employed by 
IDavis etahl ( 2013 ) and extend it to measure the SMBH 
mass using the ALMA data for NGC 1097(Project code 
= 2011.0.00108.S; PI = K. Kohno). The observations 
and data reduction are described in Section El The 
SMBH mass measurement method is explained in Sec¬ 
tion [3] Section 0] contains discussions about the eff ect o f 
dust extinction on the derived SMBH mass (Section l47Tj) . 
how the inclination angle of NGC 1097 affects the re¬ 
sult (Section l4.2j) . the dependence of SMBH mass on the 
molecular gas species used to trace the dynamics (Sec¬ 
tion [Ol). The conclusions are summarized in Section [5] 

1.1. NGC 1097 

NGC 1097 is a nearby Type- 1 Seyfert galaxy at a dis¬ 
tance of 14.5 Mpc (jTuilvl Il988h (~ 70 pc arcsecond -1 ). 
The position of the nucleus is determined by the peak 
position of the 6 cm continuum emission (Hummel et al.l 
Il987l ): iL4(J2000.0) = 02 h 46 m 18 s .96, DEC{ J2000.0) = 
—30°16 / 28 // .9. The peak position of the 860 /im contin- 
uum emission coincides with the 6 cm peak (jlzumi et all 
120131 ). Properties of NGC 1097 are summarized in Ta¬ 
ble [lj 

The SMBH mass in N GC 1097 is estim a ted t o be 
(1.2 dz 0.2) x 10 8 Mq by iLewis fe Eracleousl (120061) us¬ 
ing the empirical Mbh — & relation from [Tremaine et all 
(|2Q02h with an observed cr = 196 db 5 km s -1 . The 
uncertainty in this estimate is large, depending on the 
assumed Mbh — cr relation. The latest Mbh — cr rela¬ 
tion [log-, n (Mbh /M^) = 8.32 + 5.64 log 10 (cr/200 km s -1 ), 
[McConnell fc Mai (12013! )] would yield SMBH mass of 
(1.9 d= 0.3) x 1O 8 M 0 . Note that this relation is a fit to 
both late-type and early-type galaxies. When selecting 
only the late-type galaxies, the Mbh — & relation becomes 
log 10 (M BH /M 0 ) = 8.07 + 5.061og 10 (cr/200 km s -1 ) 


TABLE 1 

Properties of NGC 1097 


Parameter 

Value 

Reference 

Morphology 

SB(s)b 

1 

Nuclear activity 

Type 1 Seyfert 

2 

Position of nucleus 


3 

TL4(J2000.0) 

02 h 46 m 18 s .96 


DEC(J2000.0) 

—30°16 / 28 // .9 


Systemic velocity (km s -1 ) 

1253 a 

4 

Position angle (°) 

130 

1, 4 

Inclination angle (°) 

46d=5 

5 

Distance (Mpc) 

14.5 

6 

Linear scale (pc arcsec -1 ) 

70 

6 

I-band luminosity (mag) 

8.09 

7 

References. — (1) de Vaucouleurs et al. (199JLJ); (2) 


Storchi-Bergmann et al. (1993); (3) Hummel et al. (1987); 

(4) this work; (5) Ondrechen et al. (1989); (6) Tullv (1988); 
(7) ISpringob et al.l (j2007l) : 


a Systemic velocity here is a heli ocentric v e locity deter¬ 
mined with molecular lines. IKoribalski et all 12004 ) shows 
the heliocentric velocity to be 1271 km s —1 determined with 
HIPASS observation. 

([McConnell fe Mall2013l ) and the estimated SMBH mass 
becomes (1.1 db 0.3) x 1O 8 M 0 . 

The en close d mass in 40 pc radius has been studied 
by (jlzumi et alJl2013h to be 2.8 x 1O 8 M 0 , using the dy¬ 
namics fro m HCN ( J = 4 — 3) emission line. In contrast, 
iFathi et ahl (j2013h report a dynamical mass in 40pc ra¬ 
dius a s 8.0 x 1O 6 M 0 from the same data ofllzumi et all 
( 2013 ). The difference occurs because IFathi et ahl (j2013h 
assume a thin disk and extr acts the non-circular mo¬ 
tions of the gas while llzumi et~al 3 d20l3h assume a sim¬ 
ple K eplerian r otatio n. Note but the dynamical mass 
of llzumi et al ] d20l3h includes all the mass within that 
radius, not showing the intrinsic SMBH mass. A more 
detailed study of NGC 1097 is thus necessary to precisely 
measure the SMBH mass. 

2. OBSERVATIONS AND DATA REDUCTION 

NGC 1097 was observed with the band 3 receiver on 
ALMA using the two sideband dual-polarization setup 
as a cycle 0 observation (Project code = 2011.0.0108.S; 
PI = K. Kohno). The observations were conducted on 
2012 Jul 29 and 2012 Oct 19 with an hour angle from 
—4 to 2 and a total on-source time of 105.24 minutes. 
The antennas were in the Cycle 0 extended configuration 
(400m baselines) which resulted in a synthesized beam 
of 1".60 x 2".20 at a position angle —81.2° (^ 112 pc 
x 154 pc). The receivers were tuned to cover the fre¬ 
quency range from 87.275 GHz to 100.917 GHz with two 
spectral windows each in the upper sideband (USB) and 
the lower sideband (LSB). Each spectral window had a 
bandwidth of 1.875 GHz with 3840 channels. The fre¬ 
quency resolution was 0.488 MHz per channel. Obser¬ 
vational parameters are summarized in Table El The 
field of view (full width half maximum of the primary 
beam) at these frequencies was 69". The data were re¬ 
duced and imaged using CASA (Common Astronomy 
Software Applications) 4.0 with a robustness parame¬ 
ter of 0.5. We binned the data by 2 channels to im¬ 
prove the signal to noise ratio and our final resolution 
is 0.976 MHz or ^ 3 km s -1 . Molecular gas emission is 
detected over 560 km s -1 (HCN(J = 1 — 0) emission is 
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TABLE 2 

ALMA OBSERVATION PARAMETERS 


Parameter 

Date 

2012 Jul 29, Oct 19 

On-source time 

105.24 minutes 

Configuration 

extended (Cycle 0) 

Phase center: 


iM(J2000.0) 

02 h 46 m 19 s .06 

DEC( J2000.0) 

—30°16'29".7 

Primary beams 

69" 


LSB USB 

Frequency coverage (GHz) 

85.400-89.104 97.271-100.917 

Velocity resolution (km s -1 ) 

1.7 1.5 

Central frequency of 


each spectral window (GHz) 

86.338, 88.166 98.209, 99.979 


seen from 88.1524-88.3467 GHz and HCO+(J =1-0) 
from 88.7139-88.9082 GHz). The integrated intensity 
moment zero and intensity weighted velocity maps were 
made using Karma (|Goochlll996l ). These are shown in 
Figure [lj The noise in the integrated intensity maps is 
22 (mJy beam -1 km s -1 ) in the HCN(J = 1 — 0) and 
26 (mJy beam -1 km s -1 ) in the HCO + ( J =1 — 0) maps 
respectively. The peak flux is detected at 52cr in the 
HCN( J = 1 — 0) map and at 29cr in the HCO + ( J = 1 — 0) 
map respectively. The data clearly show the rotation 
dominated kinematics of the molecular gas around the 
SMBH. 

3. SUPERMASSIVE BLACK-HOLE MASS ESTIMATION 

The measurement procedure for the SMBH mass and 
its result are described in this section. We model a 
mass distribution of the galaxy with multiple Gaus- 
sians to express the combination of stellar mass pro¬ 
file and the SMBH mass. The gravitational poten¬ 
tial is derived b y fol lowing the equations described in 
iCappellari et al.l (|2002[) . which uses a Multi Gaussian 
Expansion (MGE) method. Circular velocity is cal¬ 
culated from the gravitational potential field by using 
MGE_circular_velo city c ode, which is in the JAM mod¬ 
elling packag43 of iCappellaril (j2008D . The SMBH mass 
is estimated from the total mass profile, which gives a 
velocity field best matched to the observed result. When 
comparing the derived velocity profile with the observa¬ 
tional results, we use t he KINematic Molecular Simula¬ 
tion (KinMS. iDavis et al.ll2013 ) in order to consider disk 
properties (e.g., disk thickness, position angle, and incli¬ 
nation) and the observational effect of beam-smearing. 
We show the details of each procedure in the following 
sections. 

3.1. The Mass Model 

We express the mass distribution as a summation of 
the SMBH mass and the stellar mass profile, expressed 
as the stellar luminosity distribution multiplied by a con¬ 
stant I-band mass-to-light (M/L) ratio. Assuming that 
the galaxy is axisymmetric, the stellar luminosity dis¬ 
tribution along the galaxy major axis is modeled as a 
superposition of Gaussian c omponents, using the idea of 
Multi-Gaussian Expansion (jEmsellem et al.lIl994f ). The 

6 http://purl.org/cappellari/idl 


TABLE 3 

MGE Parameters for the Stellar 
Luminosity Obtained from WFC F814W 
Surface Brightness of NGC 1097 


j Ij(Loj /pc 2 ) a j (arcsec) qj 


1 . 7451.23 0.500 0.900 

2 . 5122.72 3.50 0.700 

3 . 3725.62 1.25 0.900 

4 . 1862.81 16.0 0.700 


a See text for how we deal with the AGN and 
the starburst region. 

major axis d efined her e is at a position angle of 130° from 
Ide Vaucouleurs et al.l ( 19911 ). We use an I-band image 
observed with the Advanced Camera for Surveys Wide 
Field Channel F814W filter on Hubble Space Telescope 
(.HST ) to obtain the luminosity distribution at the major 
axis (black line in Figure [2}. We subtract the contribu¬ 
tion from the Active Galactic Nucleus (AGN) and flatten 
the starburst dominated region (both are shaded in Fig¬ 
ure [2]) in order to estimate the underlying stellar luminos¬ 
ity profile. The AGN profile is calculated by the convo¬ 
lution of a delta function with the Point Spread Function 
(PSF), measured to be 0".2 (full width half maximum, 
FWHM) from five unresolved stars in the same image. 
The PSF is also checked by using “T iny Tim” package 
(version 6.3) developed bv lKrist et al.l (|201lD . The lumi¬ 
nosity value of the delta function is determined to obtain 
the residual luminosity distribution larger than 0 at any 
radius. Note that the AGN subtraction does not gravely 
affect the SMBH mass estimation. Even the peak before 
subtraction gives a mass of less than 6.00 x 1O 6 M 0 for a 
M/L ratio 5.14, and it is at least one digit smaller than 
the SMBH mass we put in the model. The starburst ring 
region is determined to be in the range of radii from —11" 
to —7" and from 8" to 13" along the major axis. The ring 
includes younger stars than inside or outside of the ring, 
and is bluer in color. We model the stellar luminosity 
profile by calculating the least-square value with the lu¬ 
minosity distribution without the ring and the AGN. The 
model (blue line in Figure [2]), therefore, consists mainly 
of old stars, but does not include younger stars on the 
ring. See also Table [3] for MGE parameters we give for 
the data. The mass profile is simply modeled by mul¬ 
tiplying the constant I-band M/L ratio to the modeled 
stellar luminosity profile and adding the assumed SMBH 
point mass in the center. Note again that we model older 
stars, by which means we are assuming that the radial 
difference of the M/L ratio is negligible. 

3.2. Velocity Field Calculation 

With a given mass profile from Section 13.11 we calcu¬ 
late the velocity field from a gravitational po t entia l field 
derived with equations from ICappellari et al.l ((2002). We 
use MGE_circular.ve locity code, which is in the JAM 
modelling package of Cappellari (2008). The inclination 
angl e is set to be 46° (jOndrechen et al.lll989l ). See Sec¬ 
tion [472] for more details regarding the inclination angle. 

3.3. Modelling a Position-Velocity Diagram (PVD) 

We calculate a position-velocity diagram (PVD) model 
along the galaxy major axis by assuming that the ob¬ 
served molecular gas follows the velocity field obtained 
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Fig. 1. — (Left side) Integrated intensity map (moment 0) for HCN( J = 1 — 0) (upper panel, integrated for 88.1524-88.3467 GHz) 
and HCO + (J =1—0) (lower panel, 88.7139-88.9082 GHz). The rms noise level (la) in each integrated intensity map was 
22 (mJy beam -1 km s -1 ) in HCN(J = 1 — 0), and 26 (mJy beam -1 km s -1 ) in HCO + (J = 1 — 0). The map is shown in 
gray scale with contour levels 3 a to 45cr in steps of 7 a for HCN(J =1 — 0) (upper panel), and 3cr to 24<r in steps of 7 a for 
HCO + (J = 1 — 0) (lower panel). The synthesized beam size (l /x .60 x 2 /x .20 at PA = —81.2°) is shown as the ellipse at the bottom 
left of each panel. (Right side) The intensity weighted velocity map for HCN(J =1 — 0) (upper panel) and HCO + (J =1 — 0) 
(lower panel). Lower limit of the intensity is set to each map as 3cr. The velocity resolution of each image is approximately 
3.3 km s -1 . 


by our calculation. The observational effects are taken 
into account by utilizing KinMS (jDavis et al.l 20131 ). We 
convolve the model cloud distribution with our beamsize 
(1".60 x 2". 20) to express the beam-smearing, and as¬ 
sume the molecular gas disk to be an axisymmetric thin 
disk. The position angle of the galaxy major axis is set 
to be 130° (jde Vaucouleurs et al.l 1 19911 ). which is consis¬ 
tent with the kinematical position angle estimated from 
our observational data. The position angle is also con¬ 
sistent with the one calculated by Spitzer Survey of Stel¬ 
lar Structure in Galaxies (S4G. [Sheth et al.ll201Q[ ). —52°, 
and global prop erties c alcu lated via their pipeline 3 by 
iMunoz-Mateos et al.l (j2015[ ). Note that we can mostly 
avoid thjy region wit h non-circular motion pointed out 
by iFathi et all ( 2013 ) by using a PVD along the galaxy 
major axis. We assume that the streaming motion re¬ 
maining in the PVD is negligible. We can also comment 
that when considering the error propagation for a simple 
equation of v 2 / 2 = GM/r, 10 percent error in the veloc¬ 
ity could result in 20 percent error for the SMBH mass, 
which is consistent the error bar we derive from Figure [5] 


3.4. Fitting the PVD Model to the Observational Result 

The observed PVD of HCN(J = 1 — 0) is shown in 
the upper panel of Figure [3] with color filled contours. 
We set the kinematical position angle to be 130° db 5°. 
The center of the galaxy is assumed to be the peak of 
6 cm observation ({Hummel et al.l 1 198 7lh 860 /im , and 
HCN(J = 4 — 3) observation (jlzumi et al.l l2013h . We 
fit the observed PVD with PVD models calculated with 
2 free parameters, the I-band M/L ratio and the SMBH 
mass. We find the two parameters to be around M/L^5.0 
and Mbh ^ 1.0 x 1O 8 M 0 by initial robust-grid cal¬ 
culation. Then the finer grid of model parameters is 
set to be from M/L= 4.80 to 5.35 in steps of 0.01 
and Mbh = 0.5 x 1O 8 M 0 to 2.5 x 1O 8 M 0 in steps of 
0.1 x 1O 8 M 0 . We calculate an optimal rotation curve 
from the observed PVD above 3cr ^8 mJy as follows: we 
make two cuts - one in the horizontal and one in the 
vertical direction. A cut in the vertical direction gives 
us a spectrum at the pixel whereas a cut in the horizon¬ 
tal direction gives us an intensity profile at that velocity. 
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Fig. 2.— The observed I-band luminosity distribution of 
HST F814W (black line), the distribution with the AGN sub¬ 
tracted (red line), and the modeled stellar luminosity profile 
(blue line). We run the least-square value calculation at all 
radii excluding the colored regions to obtain the model pro¬ 
file. The starburst ring dominated regions (from —11 // to —7" 
and from 8 7/ to iW^ and the central area where seems to have 
the luminosity dominated with the AGN (from —0 ,x .5 to 0 /7 .5) 
are colored with gray. 


ratio and SMBH mass, by taking the parameter value 
within 99.73% confidence level (A% 2 < 9, where A% 2 = 
X 2 — Xmin )• Figure [5] show the polynomial fit to the A\ 2 
distribution for each parameter. In the left panel, the 
SMBH mass is thus estimated to be 1.40lo 32 x 1O 8 M 0 
by considering all the values below the black straight line. 
The M/L ratio is estimated to be 5.14^o‘.o4 by also con¬ 
sidering values below the black straight line in the right 
panel. Note here that the derived SMBH mass is consis¬ 
tent with th e one presumed from th e Mbh — cr relation 
reported in (Mc Connell fc Mall2Q13l) and the central ve¬ 
locity dispersion (196 ±5 km s -1 , iLewis fc Eracleousl 
l2006fF see also Section [ITT)) . 

4. DISCUSSION 

While we determine the SMBH mass and the I-band 
M/L ratio from the molecular gas dynamics, there are 
some uncertainties coming from the assumption we made 
when calculating the model. We discuss the effect of the 
observing band we us e to model the stellar luminosity 
distribution in Section 01 how the inclination angle af¬ 
fects the result (Section 14. 2D . and what if we use a differ¬ 
ent emission line to observe the molecular gas dynamics 
(Section 14.3p . 


Peak positions for both are determined with the Gaus¬ 
sian fit. We use the points when the two are consistent, 
but abandon when they do not match. In the latter case 
the spectrum is not well characterized by a Gaussian be¬ 
cause the asymmetric distribution of the molecular gas 
and the beam smearing effect is skewing the profile. 106 
points are extracted from the observed PVD (see Fig¬ 
ure n. The error bar along the velocity axis for each 
representative points of the PVD is defined to be a sum 
in quadrature of the channel width (3.284 km s -1 ) and 
the error from Gaussian fitting. After fitting a Gaussian, 
to the spectrum at each position, we determine the error 
budget to be the range of all the channels which has an 
observed value within the difference of half of the noise 
level from the maximum value of the fitted Gaussian. 

Chi-square values are calculated for each model for 
the 106 points in the observed PVD. Note that the de¬ 
gree of freedom becomes 104, because we have two free 
parameters, the SMBH mass and the I-band M/L ra¬ 
tio. Figure 0] shows the chi-square value distribution in 
the parameter space. The contour level is defined to be 
(2, 3,4, 5, 6, 8,12) x (Xmin)* The smallest chi-square value 
of 113 (reduced chi-square value is 1.09 when divided 
with the degree of freedom) is realized with parameters 
of Mbh = 1.40 x 1O 8 M 0 and the I-band M/L ratio= 5.14. 
See Figure [3] to compare three PVD models in black con¬ 
tours and lines calculated with different values of pa¬ 
rameters (Mbh = 0, I-band M/L ratio= 5.14 for the 
left column, Mbh = 1.40 x 1O 8 M 0 , I-band M/L ratio= 
5.14 for the middle column, Mbh = 7-00 x 1O 8 M 0 , I- 
band M/L ratio= 5.05 for the right column). The red 
contour in the upper panel shows the observed PVD by 
HCN( J = 1 — 0). Red dots in the middle panel represent 
the extracted points from the observed PVD. Residuals 
are plotted in the lower panels. The chi-square values are 
244, 113, and 1090 (reduced chi-square values are 2.35, 
1.09, and 10.5) for each. 

We determine the error bar for each parameter, M/L 


4.1. The Proper Stellar Luminosity Profile without the 
Dust Effect 

We estimate the expected stellar mass profile by ex¬ 
cluding the bright AGN profile and luminosity enhance¬ 
ment by the starburst ring (Section l3.ip , but we could 
not avoid the dust extinction effects, which could be im¬ 
portant for this starburst galaxy with its prominent dust 
lane around the starburst ring. One way to mitigate the 
dust extinction effects is to observe at longer wavelengths 
such as the near infrared. NICMOS on HST has a nar¬ 
row band filter F190N which observes at 1.9 microns. We 
calculate a velocity field from the assumed SMBH mass 
and the stellar mass profile derived from the luminosity 
profile of 1.9 microns, and obtain a PVD by following 
the method described in Section [3~3l We then compare 
the two PVDs calculated from the F190N luminosity pro¬ 
file and the F814W luminosity profile. Both luminosity 
profiles need to have the same PSF for a proper compar¬ 
ison. For the NICMOS data, which have a small field of 
view and no feasible stars available for a measurement 
of the PSF , we refer “Tiny Tim” package (version 6.3 
iKrist et~aDl2Qll[) and assume the shape as a Gaussian 
with its FWHM of CF.4. We convolve the I-band lumi¬ 
nosity distribution with this PSF. Figure [6] shows the 
luminosity distribution of the PSF-convolved I-band and 
1.9 microns normalized by each maximum value (upper 
panel) and the difference of two profiles (lower panel). 
This comparison already shows that the difference be¬ 
tween the two luminosity profiles is negligible. We also 
examine the S4G data at 3.6 microns to compare the stel¬ 
lar profile with the HST data but the Spitzer Infrared 
Array Camera (IRAC) PSF is too large (~ 1".8) to do 
any detailed comparison. 

We normalize the peak luminosity at 1.9 microns to 
the I-band peak, and then subtract the luminosity en¬ 
hancement of the AGN, as described in Section I3T1 We 
also ignore the starburst region and use the same M/L 
ratio to calculate the velocity field for the two stellar 
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Fig. 3. — (Upper panels) PVDs calculated with the parameter set of Mbh = 0, M/L=5.14 (left), Mbh = 1.40 X 10 8 M©, 
M/L=5.14 (middle), and Mbh = 7.00 X 1O 8 M 0 , M/L=5.05 (right). Black contours are for the modeled PVD, while color filled 
contours show the observed PVD for HCN(J = 1 — 0). (Middle panels) Each black line shows the PVD model and gray dots 
show the representative points of the observed PVD. The points are extracted from the gray contours in the upper panel by 
fitting a Gaussian to the spectrum at each position. Here we fix the lower luminosity limit to represent the observed PVD as 
3cr ~ 8.0mJy. We obtain x 2 — 113, which is the minimum value, for Mbh = 1.40 X 10 8 M© and M/L=5.14 (middle), while 
X 2 = 244 and 1090 for Mbh = 0, M/L = 5.14 (left), and Mbh = 7.00 x 1 O 8 M 0 ,M/L = 5.05 (right). We also put the reduced 
chi-square value, which is a chi-square value divided with the degree of freedom of 104, as x?ed- (Lower panels) Residual between 
the black line and plots in the above panels of each. 


mass profiles. We find that a Mbh = 1.40 x 1O 8 M 0 and 
M/L= 5.14 gives the best-fit value for an inclination of 
46°. The PVDs calculated as such are shown in Fig¬ 
ure [71 Black and gray contours in Figure [7] are the PVD 
calculated from two stellar mass profile models - these do 
not differ much between the two luminosity profiles. We 
therefore conclude that the dust extinction effect with 
the F814W filter is not too serious for the measurement 
of the SMBH mass. 

4.2. Effect from the Inclination Angle 

We discuss briefly on the difference coming from how 
we set the inclination angle. The accuracy of the incli¬ 
nation angle is critical for calculating the velocity and 
therefore crucial for the SMBH mass estimation. It is 
however not straightforward to determine the inclination 
angle when comparing observations at different field of 
view. Previous studies of NGC 1097 have determined a 
dynamical inclination a ngle of 46 zb 5° (HI obs ervations 
at 3 kiloparsec scales, lOndreche n et al.l I1989IL or 34 ° 
(Ha line profile s tudy, IStorchi-Bergman n et al.l [2003b 
iHsieh et al.l ( 2011 ) reported that the inclination angle of 
NGC 1097 is 41.7 ±0.6° using the kinematic parameters 
of 12 CO(J = 2 — 1) observed with Submillimeter Array 
( SMA ). They argue that the circumnuclear ring is nearly 
circular for the inclination of ^ 42°, by which means the 


ring has an intrinsic elliptical shape in the galactic plane, 
of which case is not symmetric to the axis. Though the 
suggested asymmetry is interesting to note, we would like 
to leave it as a further discussion, since in this work we 
assume an axisymmetric distribution for stars and molec¬ 
ular gas when calculating the circular velocity field. This 
time we assum e the galaxy inclination angle to be 46—51° 
by referring to lOndrechen et al.l (|1989h and the axis ratio 
of the HST I-band observation. We evaluate the SMBH 
mass to be x 1O 8 M 0 and the the I-band M/L 

ratio to be 5.14^o‘.o4 the inclination angle of 46°, with 
the chi-square value of 113 (1.09 when divided with the 
degree of freedom 104). We also follow the same pro¬ 
cess in Section [3] with the inclination angle of 51° and 
evaluate the SMBH mass to be 1.20^o’34 x 1O 8 M 0 and 
the I-band M/L ratio to be 5.11 ± 0.03 with the chi- 
square value of 117 (1.13 when reduced with the degree 
of freedom). See also dashed curves in Figure [5] for the 
chi-square distribution, used to determine the error bar 
for each parameter. 

We can also consider the case of the inclina¬ 
tion angle is 34 — 41° by multiplying a factor of 
sin(ii ntr insic)/ sin(46°) ^ 0.78 — 0.91 to the velocity where 
we write intrinsic as an inclination angle of the observed 
component. Under the simplified assumption that the 
SMBH mass is proportional to the square of the veloc- 
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BH mass (solar mass) 


Fig. 4.— Chi-square distribution of the grid calculation with 
two free parameters, the SMBH mass and the j-band M/L ra¬ 
tio. T he inclination angle is fixed to i — 46° (|Ondrechen et ahl 
mu. Gray dots show the parameter pair for modeling the 
PVD. The SMBH mass and the M/L ratio are fixed respec¬ 
tively with the range of Mbh = 0.50 x 10 8 M sun to 2.50 x 
10 8 M sun in step of 0.50x 10 8 M© and M/L=4.80 to 5.35 in step 
of 0.01. Contour levels are set to be (2,3,4, 5, 6,8,12) x Xmin> 
where Xmin = 113 is determined by the minimum chi-square 
value, realized with parameters of Mbh = 1.40 X 1O 8 M0 and 
M/L=5.14. We determine the error bar for each parameter 
by using Figure [5] 



Fig. 5.— Distributions of the Ay 2 = X 2 — Xmin f° r each pa¬ 
rameter, the SMBH mass (left panel) and the I-band M/L 
ratio (right panel). Solid black curve shows the distribu¬ 
tion calculated at the inclination angle of 46°. Dashed curve 
stands for the calculation for 51°. The error bar at 99.73% 
confidence level for each parameter is taken to be the values 
under solid straight line (3cr line, where Ax 2 = 9). For the 
inclination angle of 46°, the error bar for each is determined 
to be Mbh — IHOIq'J^ x 10 8 M© and M/L= 5.14 ^q'.o 4> while 
Mbh = 1.20 ±° 0 \H x 10 8 M o and M/L= 5.11 ± 0.03 for the 
inclination angle of 51°. 

ity, we can estimate the change of the SMBH mass to be 
smaller than 0.31 x 10 8 M Q , which is mostly included in 
the error bar of our result I. 40 I 0 J 2 x 10 8 M o . We there¬ 
fore consider that it is not crucial to count this error into 
our error budget. Note that, however, this galaxy could 
have a warped or a misaligned structure, which could be 




Fig. 6. — (Upper panel) The observed luminosity profile of 
I-band ( indicated with solid line, HST F814W observation) 
and 1.9 microns (dotted line, HST F190N observation). The 
luminosity value is normalized to compare the intrinsic pro¬ 
file. I-band (F814W) profile is convolved with the assumed 
poin t sp read function (PSF) of F190N. As discussed in Sec¬ 
tion ism we shade the starburst ring region and the AGN 
dominated region. (Lower panel) Residual of the two lumi¬ 
nosity distribution showing a good agreement of the two. 



Fig. 7.— Gray filled contours show the PVD calculated from 
I-band luminosity profile multiplied by the M/L ratio and 
black contours are the one from 1.9 microns luminosity profile 
multiplied by the M/L ratio. Parameters for the calculation 
is set to be Mbh = 1.40 x 10 8 M© and M/L=5.14, the best- 
fit value at the inclination 46°. PVD does not largely differ 
between the two stellar distribution obtained from different 
wavelength. 

interesting to investigate but requires a calculation for 
an asymmetric potential field. 

4.3. SMBH Mass Estimation from Other Molecular 
Species 



















































Fig. 8. — The observed PVD of HCN(J =1 — 0) (upper 
panel) and HCO + (J =1 — 0) (lower panel) is respectively 
shown with black contours. The contour level of both is from 
lcr to 4a where lcr = 4.6 mJy for HCN(J =1 — 0) and 
la = 3.2 mJy for HCO + (J = 1 — 0). The velocity structure 
of these two PVDs are in good agreement. 

Our main result is obtained using the HCN line be¬ 
cause it had the highest signal-to-noise ratio (SNR). It is 
important to measure the SMBH mass from other molec¬ 
ular species as well for consistency. We therefore repeat 
our method using the HCO + (J =1 — 0) emission line. 

We apply the fitting procedure described in Section l3~4l 
to the PVD for HCO + (J = 1 — 0), and estimate the 
SMBH mass to be (1.40 =b 0.30) x 10 8 M© and the I-band 
M/L ratio to be 5.15 db 0.03 with a galaxy inclination of 
46°. These derived values are consistent with the mea¬ 
surement using HCN(J = 1 — 0). From FigureEl we see 
that the observed PVDs of two molecular gases are in 
good agreement, indicating that the fitting parameters 
will be consistent between the two. 

Reaching the velocity structure from multiple molecu¬ 
lar species is one of the particular benefit of the SMBH 
mass measurement with millimeter/submillimeter wave¬ 
length observation, which enable one to observe more 
than two molecular species at the same frequency band. 

5. CONCLUSION 

We derive the SMBH mass in NGC 1097 to be 
1.40l8;g x 1O 8 M 0 by using dense molecular gas dy¬ 
namics traced with HCN(J = 1 — 0) and HCO + (J = 
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